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& ■ Abstract 
O 

The susceptibility of an overdamped Markov system fluctuating in a bistable potential 
• i— i . 

of general form is obtained by analytic solution of the Fokker-Planck equation (FPE) for 
low noise intensities. The results are discussed in the context of the LRT theory of 
stochastic resonance. They go over into recent results (Gang Hu et al Phys. Lett. A 
172, 21, 1992) obtained from the FPE for the case of a symmetrical potential, and they 
coincide with the LRT results (Dykman et al, Phys. Rev. Lett. 65, 2606, 1990; JETP 
Lett 52, 144, 1990; Phys. Rev. Lett. 68, 2985, 1992) obtained for the general case of 
bistable systems. 
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The phenomenon of stochastic resonance (SR), where a periodic signal in a system 
and the signal-to-noise ratio are enhanced by external noise, has attracted much attention 
recently and by now has been observed in several physical systems (see [1] and references 
therein). The occurrence of SR is related to the nonlinearity of a system. At the same 
time, for small enough amplitude A of the periodic force causing a signal, the amplitude 
of the signal is proportional to A, and therefore many aspects of SR can conveniently and 
effectively be analysed [2-6] in terms of standard linear response theory (LRT) [7]. The 
LRT analysis amounts to finding the susceptibility of the system x(f2) an d its dependence 
on the noise intensity D and on the frequency Q of the force. The susceptibility is defined 
in a standard way [7] as the coefficient of proportionality between the periodic term in the 
ensemble-averaged value of the coordinate 5{q(t)) (the signal) and the force Aexp(—iQ,t), 

6{q(t))=x(n)Aexp{-int) (1) 

(in the case of a nonstationary system the definition (1) has to be modified). The am- 
plitude of the signal is equal to |x(fi)|A, and SR arises when increases with D in 
some interval of D (and then decreases again for higher D, so that the dependence on 
noise intensity is of the form of a resonant curve, thus justifying the term "stochastic 
resonance" ) . 

The LRT formulation is not confined to particular systems or particular types of noise 
and, indeed, was used to predict and observe SR [4, 6] in systems drastically different 
from that performing Brownian motion in a symmetric double-well potential which is 
usually considered as the prototype of SR-displaying systems [1]. One of the advantageous 
features of this formulation is that, in the case of systems at thermal equilibrium or quasi- 
equilibrium (in which case it is temperature T that corresponds to the noise intensity D), 
the susceptibility x(fy i s related, via the fluctuation-dissipation relations, to the spectral 
density of fluctuations (SDF) of the coordinate q(t) [7]. Therefore, if the SDF is known 
from the experimental data, it is straightforward to predict whether the system will or 
will not display SR when driven by a force of a given frequency even without establishing 
the underlying physical mechanism; we note that the corresponding analysis does not 
assume that the system is Markovian or whatever - the only thing assumed is that it is 
in thermal equilibrium. 
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The theoretical evaluation of the susceptibility for some simple model noise-driven 
systems that display SR, and for systems performing Brownian motion in a double-well 
potential in particular, has been done in Refs. [2 - 6] . The analysis of bistable systems 
could be done for small noise intensities where there are two strongly different time scales 
that characterize the motion of a system: the relaxation time t r over which the system 
approaches a stable state (the longest of these times), and the reciprocal probabilities 
W~m of fluctuational transitions n — > m (n,m — 1,2) between the states, 

Wt r < 1, W = W 12 + W 21 (2) 

Two approaches have been used. The first of these [2, 3] is based on the idea that, for low 
noise intensities, a system spends most of its time localized close to the stable states, and 
only occasionally do there occur fluctuational transitions between the states; therefore the 
susceptibility is a sum of the partial susceptibilities Xn(^) from small- amplitude vibrations 
about the stable positions q n (n = 1,2), weighted with the populations w n of the stable 
states, and of the susceptibility related to the fluctuational transitions, Xtr(fy, which 
describes the effect of the modulation of the transition probabilities and hence of the 
populations of the stable states by a periodic force (cf. [8]): 

X(P.) = Yl WnXniti) + XtrM, W n = W mn /W (3) 
n=l,2 

The evaluation of Xn reduces to the solution of a linearized problem, and the explicit 
expression for Xtr has been obtained for systems driven by Gaussian noise in the range 
Vlt r <C 1 (xtr oc W/VL for W <C fi, cf. Eq.(18) below, and it becomes very small for 
VLt r . ~ 1 in the range (2)). 

A different approach has been suggested in Ref.5 for the analysis of overdamped Brow- 
nian motion in a symmetrical bistable potential U(q), 

q + U'(q) = f(t), (f(t)f(t>)) = 2D5(t-t>) (4) 
The basic idea of this approach is to reduce the Fokker-Planck equation 

dP(q,t) d(U'(q)P(q,t)) ^ d 2 P(q,t) 

= dq + D ^q^~ (5) 
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that corresponds to the Langevin equation (4) to an eigenvalue problem and to express 
the susceptibility in terms of the matrix elements of the coordinate q on the corresponding 
eigenfunctions. It was then possible to evaluate the susceptibility explicitly in the range of 
noise intensities D where (2) holds. A systematic perturbation theory in the amplitude of 
a sinusoidal force has been developed for systems described by the Fokker-Planck equation 
and applied to SR in an effectively two-state approximation in Ref [9]. 

The purpose of the present paper is to demonstrate that the two approaches [2, 3] 
and [5] give identical results and to generalize the approach [5] to the case of an arbitrary 
double- well potential U(q). 

The reduction of Eq.(5) to a Hermitian eigenvalue problem is standard [10, 11]: one 
seeks the solution in the form 

oo 

P(q,t) = Z- 1 / 2 e W (-U(q)/2D)J2\j)eM-Hm, \j)=^U\ </)> (6) 

3=0 

/oo 
exp(-U(q)/D)dq 
-oo 

and then arrives at a Schrodinger-type equation for the eigenfunctions q): 

- D 2 r(j; q) + V(q)W, q) = DX{j)m ?), V(q) = \ (U' (q) f - l -DXJ" \q) (7) 

The functions q) are orthonormal, (j\f) = Sjj> (the form of the expansion (6) is a bit 
different from that used in [5]). Following the arguments given in Ref. 5 one can arrive 
at the following expression for the susceptibility of the system (4) with respect to a force 
Aexp(— iQt): 

x(n) = -2£(A(j) - in) _1 (o|?U>0'l^|o> (0'W> = JZ dqm q)L ^ 3 ' ] q) ) (8) 

The eigenvalues X(j) and the eigenfunctions ip(j; q) for low noise intensities D were anal- 
ysed in Ref. [11]. The qualitative results can be easily understood by noticing that, for 
a bistable initial potential U(q) and small D, the potential V(q) that appears in Eq. (7) 
has 3 minima of nearly the same depth placed close to the stable positions q ij2 and an 
unstable stationary point q s (local maximum of U(q); U'(qi : 2) = U'(q s ) = 0, t/"(<?i,2) > 
0, U"(q s ) < 0). Close to the minima V(q) is parabolic: 
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V(q)^--DU^ + -(U^) 2 (q-q n ) 2 , \q - q n \ « \q n - q s \, U^ = U"(q n ) (9) 

V(q) » ±D\U'J\ + \(U' s ') 2 (q - q s )\ \q - q s \ « |g 1>2 - q s \, U'J = U"(q s ) 

The lowest "energy levels" A n (0) of the motion within the wells n = 1, 2 of the potential 
V (the overbar is used to indicate intrawell states , and the subscript enumerates the wells 
at the corresponding stable positions q n ) are degenerate, with Ai^O) = to the lowest 
order in D, giving rise to tunnelling splitting of these levels [11]. The low-lying excited 
intrawell states \j) n = ip n (j; q) are those of quantum harmonic oscillators, and 

K{k) » kU" (10) 

(and of course, high excited states cannot be prescribed to a particular well). 

The exact wave function of the ground state (that allows for the tunnelling splitting) 
is known to be 

^(0; q) = Z- 1 / 2 exp(-U{q)/2D), A(0) = (11) 

It has sharp maxima at qi :2 (of width 5q ~ (D /U"{q n )) 1 ^ 2 ) and a minimum at q s . Near 
the maxima, the function (11) is close to the intrawell ground-state oscillatory functions 
V^(0; q) weighted with the factors Z- 1 / 2 (2ttD/U^) 1/4 x exp(-U(q n )/2D). This makes it 
possible to evaluate the contribution to the susceptibility (8) from the excited intrawell 
states, XintraiS^)- Since, for a harmonic oscillator, the only nonzero matrix elements of the 
operator d/dq are those between adjacent states, only the first excited intrawell states, 
V'O'i q) ~ V'nClj l) ( n — 1)2) will contribute to XintraiS^)- Allowing for the explicit form 
of the matrix elements, known from quantum mechanics [12], and for Eqs.(9), (10) we 
get, with an accuracy to small corrections ~ D: 

XintraW » £ W n {U' r [ - 1^ (12) 
n=l,2 

Here, w n is the population of the nth stable state, 
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W r , 



exp(-AU 3 _ n /D)(Ul n ) 



II \l/2 



r , AU m = U{q s ) - U(q r , 



(13) 



E m=1)2 exp(-AC/ m / J D)(C/-)V2^ 
and we have evaluated the integral (6) for Z by the steepest descent method. The overlap 
integral of the function ip(0; q) with (non-oscillatory) wave functions of highly excited 
states, as well as with the wave functions that correspond to the low-lying states in 
the well of the potential V(q) centered at q s , are exponentially small, and therefore the 
contribution from the corresponding matrix elements to the susceptibility (8) can be 
ignored. 

The only other important contribution to x(^) besides (12) is that from the matrix 
elements calculated between the function (11) and the function ^(1; q) corresponding to 
the excited state that results from the tunnel splitting of the degenerate intrawell state, 
whence 



^(fi) ~ Xintra{ty + Xtun(ty, Xtun{ty 



A(l) |<0|?|1>| S 



(14) 



D A(l)-iQ 

Here, we have implemented the interrelation [11], well-known from quantum mechanics, 
between the matrix elements of the operators d/dq and q in terms of the energy difference 
between the states. The problem of evaluating xtun comes to the evaluation of the matrix 
elements of the coordinate. The latter can be done by making use of the explicit form of 
the wave functions |0) (11) and |1). The function |1) as well as the value of A(l) were found 
in [11]. An alternative way to find them is based on applying to the present problem the 
approach suggested in [12] to calculate tunnelling level splitting in a symmetrical double- 
well potential. It is easy to see that, in the spirit of [12], the wave function ip(l; q) can 
be sought in the form 



^(1; q)=A ^(0; q) + q) 



(15) 



m q) 



-^(0; q), 
m q)(2\U' s '\/nDy/ 2 g dq> exp 
^(0; q), 



U{q')-U(qs) 
D 



? S -?»(» 1/2 

\q-q s \Z{D/\U>:\yt 
q-q s ^(D/\U'J\y/* 



The coefficients A ^ can be obtained from the conditions of normalization of ^(1; q) and of 
orthogonality of ip(l; q) and ■*/>((]; q), whereas the eigenvalue Ai is given by the expression 



A(l) = -D^O; q s ) (^^j [f^dqm 9^(1; q) 



(16) 



that immediately follows from (7), (11). 

It is seen from Eqs. (8), (14) - (16) that the susceptibility XtuniS^) due to the transitions 
to the tunnel-split state — > 1 is of the following form 



XtunM = w lW2 (q 2 - qi) 2 ^- (A(l) - in)' 1 (17) 



A(l)=vr- 1 £ 

71=1,2 



UZ\UZ\exp(-AU n /D) 



In the particular case of a symmetrical potential U(q) Eqs.(l), (12), (14), (16) go over 
into the result obtained in [5]. 

It is interesting to compare the result of the present calculations with the expression 
(3) applied to the model (4). The explicit form of the susceptibility for this model was 
discussed in [2, 3]. The intrawell susceptibilities Xn(^) = (U„ — if2) -1 follow from Eq.(4) 
linearized about the stable positions q± j2 . The expression for the susceptibility due to 
interwell transitions is of the form: 

Xtr(fi) = w lW2 (q 2 - qi ) 2 ^ (W - til)- 1 (18) 

where w n is the population of the nth state, and W is the sum of the probabilities of 
the transitions between the states (cf. (2)). Allowing for the explicit expressions for the 
transition probabilities, one finds that the intrawell contribution Xintra(ty obtained above 
is precisely equal to that given by (3), and that A(l) = W (cf. [10, 11]) and therefore 
Xtun{0) — Xtr{ty, and the susceptibilities as a whole coincide with each other. So they 
do also in the particular case of a symmetric double- well potential discussed in [5]. The 
explicit expression for the spectral density of fluctuations considered in [5] coincides with 
that discussed in [2,3] for the model (3) as well. We notice that, in the general case of an 
asymmetric potential, the intensity of the transition- (or, equivalently) tunnelling-induced 
peak in the susceptibility (17), (18) is proportional to W\W 2 oc exp(— \U(qi) — U(q 2 )\/D), 
and it drops down extremely sharply with the difference in the depths of the wells. 

It follows from the above discussion that two, seemingly different, approaches to the 
analysis of the susceptibility of a double-well system, one based on a simple physical 
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picture of the motion and the other based on the solution of the Fokker-Planck equation 
when the latter applies, give identical results. This can be regarded as an extra indication 
of the usefulness of linear response theory in the context of stochastic resonance. 
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